Flux pinning in (1111) iron-pnictide superconducting crystals 
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Local magnetic measurements are used to quantitatively characterize heterogeneity and flux line 
pinning in PrFeAsOi- y and NdFeAs(0,F) superconducting single crystals. In spite of spatial fluctu- 
ations of the critical current density on the macroscopic scale, it is shown that the major contribution 
comes from collective pinning of vortex lines by microscopic defects by the mean-free path fluctu- 
ation mechanism. The defect density extracted from experiment corresponds to the dopant atom 
density, which means that dopant atoms play an important role both in vortex pinning and in quasi- 
particle scattering. In the studied underdoped PrFeAsOi_ y and NdFeAs(0,F) crystals, there is a 
background of strong pinning, which we attribute to spatial variations of the dopant atom density on 
the scale of a few dozen to one hundred nm. These variations do not go beyond 5 % - we therefore 
do not find any evidence for coexistence of the superconducting and the antiferromagnetic phase. 
The critical current density in sub-T fields is characterized by the presence of a peak effect, the 
location of which in the (B, T)-plane is consistent with an order-disorder transition of the vortex 
lattice. 

PACS numbers: 74.25.Sv; 74.25.Uv; 74.70.Xa ; 74.25.Wx ; 74.62. En 



I. INTRODUCTION 

The characterization of the physical properties of new 
superconducting materials such as the recently discov- 
ered iron-pncitide superconductors 1 - - — requires a good 
knowledge of sample morphology and microstructure. 
The measurement and interpretation of thermodynamic 
quantities such as the magnetization, the magnetic 
torque*^ or the specific heat, or transport properties such 
as the resistance or irreversible magnetization, may be 
complicated by material inhomogeneity on mesoscopic 
or macroscopic length scales. On the other hand, mi- 
croscopic disorder is well-known to be beneficial for vor- 
tex line pinning and high critical currents. Finally, from 
the defect-vortex interaction, one might hope to extract 
information on electronic scattering mechanisms in the 
iron-pnictide superconductors, as well as on the premise 
of phase co-existence. In underdoped pnictides especially, 
it has been argued that the coexistence of the low-doping 
anti-ferromagnetic state and the superconducting state 
at higher doping levels may affect physical properties^ 

Vortex pinning and the critical current density in the 
iron pnictide superconductors has mainly focussed on 
the so-called "122" compounds, since large single crys- 
tals of these are available. Most notably, magnetic flux 
penetration in Ba(Feo. 93000.07)2^82 has been studied us- 
ing magneto-optical imaging by Prozorov et aZ. i n i 12 The 



same authors reported on the irreversible magnetization 
and flux creep in this compound, and found qualita- 
tive agreement with collective creep in the so-called bun- 
dle regime^ 3 - The non-monotonous behavior of the sus- 
tainable current as function of magnetic field was inter- 
preted in terms of a crossover to plastic creep A sim- 
ilar behavior was found for crystals with different dop- 
ing levels j 12 ' 14 the overall behavior of the critical cur- 
rent density as function of doping was attributed to the 
changing density of structural domain walls, that act 
as strong pinning centers^ Yamamoto et al. obtained 
similar results on the same Ba(Feo.gCoo.i)2As2, but at- 
tributed the temperature- and field-dependent features 
of the critical current density to an inhomogeneous disti- 
bution of Co atoms ^ Very large critical currents, as well 
as a non-monotonous width of the irreversible magneti- 
zation loops correponding to a peak-effect in the criti- 
cal current 1 ^— were measured by Yang et al. in sin- 
gle crystalline Bao.6Ko.4Fe2As2r^ who concluded to the 
presence of small-sized normal state regions in their sam- 
ples. Finally, irreversible magnetization and flux creep 
measurements were conducted on SmFeAsO0.9F0.1- and 
polycrystalline NdFeAsOo.82Feo.i8 ) 24l2S members of the 
"1111" family of compounds. In all the above cases, the 
critical current at low fields was characterized by a peak 
and negligible magnetic relaxation, followed by more pro- 
nounced thermally activated flux motion at higher fields, 




which was found to be in qualitative agreement with the 
collective creep theory*^ However, no quantitative anal- 
ysis of the data has been performed, and no definite con- 
sensus as to the defects at the origin of flux pinning has 
been established. 

The aim of the present paper is the identification of 
defects responsible for flux pinning in single crystals of 
the (Re)FeAsO "1111" family of superconducting com- 
pounds. The microstructure is characterized by the un- 
dulation of the FeAs layers and the presence of sparse 
nanometer-sized defects, both of which do not seem to 
influence flux pinning. The largest contribution to the 
critical current j c is shown to arise from the dopant 
atoms, which act as scatterers for quasi-particles in the 
vortex cores. One therefore deals with pinning by lo- 
cal variations of the mean- free path {5k mechanism). 
The temperature- and field dependence of j c is very well 
described by collective flux pinning in the single-vortex 
limit, but superposed on a strong pinning contribution 
arising from small fluctuations of the doping level on the 
scale of dozens of nm. 



II. EXPERIMENTAL DETAILS 

PrFeAsOi-j, crystals (with the P4/mmm structure) 
were grown at 1300°C and 2 GPa from pressed pel- 
lets consisting of the starting materials PrAs, Fe, and 
Fe2C>3, in the nominal composition PrFeAsOo.6' 2Sl27 The 
typical size of the crystals is 100 x 100 x 30 /im 3 ; the 
average final composition corresponds to y ~ 0.1. A 
number of monolithic crystals from this batch has been 
previously used for the measurement of the superfluid 
density^ the field of first flux penetration^ and the 
electrical resistivity in the vicinity of the upper critical 
field B c2 = $o/2tt£ 2 (with £ the coherence length) ,22 The 
superconducting properties of the compound are there- 
fore completely characterized. The temperature depen- 
dence of the in-plane penetration depth X a b(T) (for cur- 
rents parallel to the ab plane) is well described by a simple 
two-gap model, without any nodes of the order param- 
eter. The magnitude of the low temperature penetra- 
tion depth is Xab(0) = 280 nm^l Table [Q gathers the 
superconducting parameters of PrFeAsOi_j,, including 
the characteristic energy e = ^o/^^MoA^ (correspond- 
ing to 4£ 2 times the condensation energy), the low-field 
anisotropy ratio s\ = A a b/A c , and the depairing current 
jo = 4e /v / 3$o£ (jm = 4tt x 10" 7 Hm" 1 ). 

Several PrFeAsOi-y crystals were prepared for Trans- 
mission Electron Microscopy (TEM). In each case, two 
crystals, of lateral dimensions ~ 100 /im, were glued be- 
tween 0.5 mm thick Si platelets; these were then thinned 
down until the crystals were flush with the edges. Further 
thinning yielded sections parallel to the c-axis, suitable 
for TEM. FigureQJi, a high resolution image of one of the 
sections, shows clear contrast corresponding to the FeAs 
planes, with some undulation. The presence of 5 - 10 nm 
sized defects, possibly secondary phase precipitates, is 




FIG. 1: Transmission electron microscopy images of single 
crystalline PrFeAsOi-y (a) High resolution bright field image, 
revealing the undulation of the FeAs layers, (b) Bright field 
micrograph of a zone containing a nm-size inclusion (indicated 
by the white arrow), (c) Bright field image revealing contrast 
due to a line dislocation (indicated by the arrow). 



also observed (Fig. QJ>) . These defects are separated by a 
distance of the order of several dozen to several hundred 
nm, depending on location. Finally, Fig. [T] shows con- 
trast associated with the presence of a linear dislocation 
core, occasional examples of which were found. 

Three PrFeAsOi_ y crystals (# 1, 3, and 7) were char- 
acterized by X-ray diffraction usin g 28.3 keV (0.43811 A) 
radiation on the CRISTAL beam at the SOLEIL syn- 
chrotron. Images of diffraction spots were collected on a 
2D CCD detector when the sample was rotated around 
an axis in the ab plane. From the 360 measured images, 
successively collected during 1 s after a progression of 
1° of the crystal rotation, layers of reciprocal space were 
numerically reconstructed. Three such sections, contain- 
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FIG. 2: (color online) Cuts through reciprocal space of the 
PrFeAsOi-y compound, reconstructed from synchrotron ra- 
diation X-ray diffraction on crystal # 7 (the same crystal as 
in Ref.—). (a) the [hkO] plane (b) the [Okl] plane; (c) the [hOl] 
plane. 



ing the origin, are shown in Fig. O The [hkO] section 
reveals very good translational order in the basal plane. 
However, the fulfilment of the Laue condition over ex- 
tended streaks in the [hOl] and [Okl] planes shows that 
crystalline order along the c-axis is not as good. The pro- 
nounced elongation of both low- and higher order nodes 
in the [001] direction indicates that this disorder more 
than likely originates from the undulation of the planes 
observed in TEM. Other kinds of c-axis disorder such 
as stacking faults or anti-phase boundaries would have 
yielded a larger broadening of nodes outside the [hkO] 
plane, as compared to the lower order nodes. From the 
elongation of the nodes at [100] and [010], we estimate 
the buckling of the layers to result in a variation of their 
orientation of up to 5°. The same results were obtained 
for all studied crystals. 

The NdFeAsOn.flFn i crystals used in this study are the 
same as that of Ref. [29|; they were synthesized at high 
pressure in a cubic, multi-anvil apparatus. The crystals, 
extracted from a polycrystalline batch, had dimensions 
210 x 320 x 30 ^m 3 (# 1) and 150 x 200 x 50 ^m 3 (# 
2), and critical temperatures T c = 34.5 ± 1.5 K 1) 
and 37.5 ± 1 K (# 2). The superconducting parameters 
of NdFeAs(Q,F) of this particular doping level have been 
studied in Refs.[29andl3Cl and are summarized in Tablefl] 

In order to obtain the value and local distribution of T c 
and j c , flux penetration into the superconducting crys- 
tals was imaged using the direct magneto-optical imag- 
ing (MOI) method^ Crystalline inhomogeneity in the 
vicinity of the critical temperature was characterized us- 
ing the Differential Magneto-Optical (DMO) method^ 
In MOI, a ferrimagnetic garnet indicator with in-plane 
anisotropy is placed on top of the sample under study, 
and observed using a polarized light microscope. The 
presence of a non-zero perpendicular component B j_ of 
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FIG. 3: (color online) (a) DMO images of the screening of an 
ac field of 1 Oe by PrFeAsOi- H (bi)crystal # 1, at various 
temperatures close to T c . Screening first appears in the upper 
left-hand corner by the crystal (indicated by the arrow). The 
crystal progressively transits to the superconducting state be- 
tween 38 and 32 K. (b) Local permeability, defined as the ratio 
(J(r, T) - I(r, T < T c )/(J(r, T » T c - J(r, T < T c ) of the 
local luminous intensities I(r,T), in the zones 1-4 indicated 
in (a). 



the magnetic induction is revealed, by virtue of the Fara- 
day effect of the garnet, as a non-zero intensity of re- 
flected light when the polarizers of the microscope are 
(nearly) crossed. Thus, light areas in the MO images cor- 
respond to areas of high perpendicular induction, while 
dark regions have small or zero B±. In DMO, magneto- 
optical images taken at applied fields H a and H a + AH a 
(with AH a = 1 Oe) are subtracted; the procedure is re- 
peated 100 times, and the subtracted images averaged. 

The local critical current density of the investigated 
crystals was obtained by calibrating the luminous inten- 
sity of the MOI images, so as to obtain a map of the 
local induction. j c was then determined as twice the gra- 
dient of the local flux density, measured over an interval 
of length 20 ^m perpendicular to the sample boundary, 
and averaged over a width of 20 /im, (parallel to the sam- 
ple boundary). This procedure is justified in that, given 
our crystals' aspect ratio, flux profiles at the sample sur- 
face are nearly linear^ In what follows, the area over 
which j c was measured was chosen such that (B) ss 300 
Oe over the 20 x 20 ism 2 region. 

Further measurements were carried out using a micron- 
sized Hall probe array, tailored in pseudomorphic 
GaAlAs/GaAs heterostructure^i Tne 10 Hall sensors 
of the array had an active area of 3 x 3 /im 2 , while an 
eleventh sensor was used for the measurement of the ap- 
plied field. The Hall probe magnetometry technique is 
complementary to magneto-optical imaging in that it has 
greater sensitivity and can be used up to substantially 
higher magnetic fields; on the other hand, it only allows 
the measurement of B± along the array of sensors, and 
not over the entire two-dimensional sample surface. 



compound 


A„ 6 (0) 


£a(0) 




(0) 




so(0) 




io(o) 




Gi 


PrFeAsOo.sr^ 


280 nm 


0.4 


1.8 


nm 


3.2 x 


w- VJ 


Jm" 1 


2 x 10 12 Am" 2 


3 


x 10 _,i 


NdFeAsOo.gFo.i— ^ 


270 ± 40 nm 


0.25 


2.4 


nm 


3.5 x 


lO" 12 


JnT 1 


1.6 X 10 12 Am" 2 


3 


x 10" 3 



TABLE I: Superconducting properties of the crystals used in 
this study. 



III. RESULTS 
A. PrFeAsOi y 

Spatial inhomogeneity of the critical temperature in 
single crystals was investigated by the DMO images near 
the transition. A typical example is shown in Fig.[5Ji, de- 
picting four DMO images, acquired with a AH a = 1 Oc 
modulation in the absence of a static field, at various tem- 
peratures spanning the normal-to-superconducting tran- 
sition. In this particular case, diamagetic screening first 
appears at T ~ 38 K in the upper left-hand corner of 
the crystal. Magnetic flux is progressively excluded from 
the crystal bulk, until the largest part is fully screened 
at T = 34 K. However, the small grain at the bottom 
is only fully screening at T = 31 K. Fig. [5}} shows the 
ac permeability, determined from the luminous intensi- 
ties J(r,T) as T MO = (/(r,T) - J(r,T « T c )/(/(r,T » 
T c — I(r,T <C T c ), for four regions indicated in the last 
panel of Fig. [3]}. It is seen that, locally, the crystal shows 
sharp transitions to the superconducting state. However, 
a global measurement (e.g. by a commercial magnetome- 
ter) would clearly result in a broadened transition. 

Local values of the critical current density j c are ob- 
tained from the MO imaging of the largest grains in 
polycrystalline conglomerates, or from the flux distribu- 
tion in monolithic crystals, such as depicted in 0^, for 
PrFeAsOi_y crystal # 7. The magnetic flux distribu- 
tions in such crystals are characteristic of the Bean crit- 
ical state | 33 i 35 ~ — Fig. [5] shows an example of profiles ob- 
tained across the central part of crystal #7atT=llK. 
Due to the relatively large thickness-to- width ratio of the 
crystal, d/w ~ 0.3, flux profiles resemble straight lines; 
j c w T;dB±/dx can be straightforwardly obtained from 
the flux density gradient^ 3 - 

Resulting values of the critical current density in four 
areas of PrFeAsOi_j, crystal # 7 are shown in Fig. as 
function of temperature. The inset to the Figure reveals 
the inhomogeneity of T c for this particular crystal; the 
regions in which j c was measured are also indicated. It is 
found that j c = 3± 1 x 10 9 Am~ 2 at the lowest measured 
temperature. The temperature dependence j c {T) de- 
pends on location. Low j c areas show a smooth decrease 
with temperature, whereas regions where j c is higher fea- 
ture a crossover in the temperature dependence. Similar 
behavior is found in all investigated PrFeAsOi_j, crys- 
tals, see Fig. [5J We shall, in section IIV1 attribute this 
behavior to the additive effect of weak collective pinning 
by oxygen dopant atoms, yielding a strong temperature 
dependence, and strong pinning, with a weak temper- 
ature dependence, coming from disorder of the doping 
level on the scale of 10 - 100 ran . 

Measurements in higher magnetic fields were per- 
formed using the Hall array magnetometry technique. 
Typical results for the self-field, defined as H s = B±/^io — 
H ai measured on the central part of the top surface of 
crystal =^ 7, are shown in Fig. [7] The screening cur- 
rent density is proportional to the difference AH S mea- 
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FIG. 4: (color online) (a) Direct MOI of the screening of an 
dc field by PrFeAsOi- y crystal # 7 (the same crystal as in 
Ref.— ), at 7.1 K. Shown are a polarized light image of the 
crystal; and MOI for H a = 300 Oe, 500 Oe, and the trapped 
flux in zero field (after the application of 500 Oe). (b) Flux 
profiles, measured from top to bottom across the central part 
of this crystal, for increasing values of the applied magnetic 
field H a , and a temperature of 11 K. 



sured on the decreasing- and increasing field branches, 
respectively. A clearly non-monotonous field-dependence 
of the critical current is observed, with the sustainable 
current density j rapidly decreasing as the H a is first in- 
creased, followed by an intermediate regime of constant 
j. Fig.fTlTa') shows that the low-field behavior, a plateau 
up to B* , followed by a power-law decrease ~ B~ 5 / 8 , is 
archetypal for a strong pinning contribution to the crit- 
ical current. However, at intermediate fields, around 0.1 
T in Fig. [TTT a). j c does not vanish, but saturates at a 
value jf? v ~ 2 - 3 x 10 9 Am -2 at low temperature. The 
temperature dependence of the zero-field- and interme- 
diate (constant) values of the critical current are plotted 
in Fig. [HI One sees that the contribution is spatially 
rather more homogeneous, and also that it corresponds 
to the critical current measured in the most weakly pin- 
ning areas of the crystals. Below, we shall attribute this 
contribution to weak collective pinning by dopant atoms. 
The strong pinning contribution j c (0) strongly depends 
on the location at which it is measured, and it is respon- 
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FIG. 5: (color online) Critical current density j c of 
PrFeAsOi_ B crystal 7 at H a = 300 Oe, as function of temper- 
ature. j c was determined from the slope of the local magnetic 
flux density, measured in the different regions of the crystal 
indicated in the Inset. The lower drawn line is a fit to Eq. 
with rid = 1.5 x 10 27 m~ 3 , the upper drawn line shows that 
the low temperature j c (T) in the strongly pinning regions is 
the same. Inset: Polarized light image of the crystal, and five 
DMO images of screening of AH a = 1 Oe at T = 34 - 30 
K (in steps of —1 K). The regions 1-4 over which j c was 
determined are indicated by drawn squares. 



sible for the larger measured critical current densities. 

Finally, we turn to higher applied magnetic fields. It 
is observed that the hysteresis loops open up at a field 
B on , corresponding to the increase of j at the so-called 
"fishtail" or peak-effect^M^^-. 2 ^ The B on (T) data 
are collected, together with the irreversibility fields de- 
termined from the appearance of a third harmonic com- 
ponent in the ac-response^, in Fig. [12] 

As in previous studies on other iron pnictide 
superconductors ) 11 ! 14 ! 22 ! 23 ! 25 the local flux density in 
Tesla fields is observed to decay with time, with a typi- 
cal relaxation rate S = (din B±/ dint) ~ —0.05 for fields 
below H on and S ~ —0.03 for H a > H on . As in other 
studies^ magnetic relaxation was not observed to affect 
the low-field MO data, it therefore does not affect the 
measured temperature dependence of the critical current 
density in what follows. 



B. NdFeAsOi.^F^ 

Fig. [S^t shows magneto-optical images of flux penetra- 
tion into NdFeAsOo.gFo.i crystal # l. The sample turns 
out to be a bicrystal, with a similar spread in T c as ob- 
served in PrFeAsOi_j,. As shown by Fig. [5)3, flux dis- 
tributions inside the crystalline grains are well-described 
by the Bean critical state models Local values of the 
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FIG. 6: (color online) PrFeAsOi_ H : critical current versus 
temperature, for two different crystals. Open and closed cir- 
cles represent the upper and lower bounds of the j c of crystal 
#7,— the closed squares represent the zero-field j c of crys- 
tal #3, while the open squares depict the magnitude of the 
critical current in the field regime between B* , above which 
the strong pinning contribution becomes negligible, but below 
the onset field of the "fishtail" effect, B on . This value of j c 
is attributed to weak collective pinning by oxygen vacancies 
in the single vortex limit; the drawn line represents a fit to 
Eq. (O, with n d = 1.5 x 10 27 m~ 3 . The dashed line is a fit 
to Eq. (|f Op summing strong and weak pinning contributions 
to the critical current at zero field, with parameter values 
n d = 1.5 x 10 27 m -3 , m = 1 x 10 21 m -3 , and f PlS = 0.1e 



critical current density at B± sa 300 Oe were obtained 
in the same manner as described above. Results for the 
three regions outlined in the center panel of Fig. [5ji are 
rendered as function of temperature in Fig. [5] together 
with results obtained by Hall probe magnetometry over 
the central regions of crystals #1 and 2. Field-dependent 
results are shown in fig. Illfb ). The overall behavior re- 
calls that reported in Ref. [25|, and is very similar to that 
observed in PrFeAsOi_ y : high critical current areas cor- 
respond to a large local contribution of strong pinning, 
whereas the lower v measured at intermediate fields 
much above B* corresponds to the critical current den- 
sity in the more weakly pinning areas of the crystals. 
In contrast to PrFeAsOi_j,, the strong pinning contri- 
bution outweighs jf y by a factor 2-3. NdFeAsO(0,F) 
crystal # 2 shows a clear "fishtail" or peak-effect, the 
corresponding B on (T) values are plotted in Fig. [T3J A 
hint of a peak-effect is also observed in crystal # 1, but 
the relative increase of the sustainable current density is 
much more modest than in the other inve stig ated sam- 
ples, with data resembling those of Ref. [25l Finally, 
Fig. 1 121 shows that the irreversibility field measured from 
the onset of screening^ coincides with that determined 



(i 



from the onset of a third harmonic response in ac Hall- 
probe array magnetometry. Moreover, the irreversibility 
field B„ r (T) for NdFeAs(0,F) and PrFeAsOi_ y crystals 
with the same T c are, within experimental accuracy, iden- 
tical. 



IV. DISCUSSION 



A. Weak collective pinning 



We start by analyzing the critical current contribution 
in terms of the weak collective pinning theor y 13 i 39 i 40 
The vortex lattice order is characterized by the transverse 
and longitudinal displacement correlation lengths 



(\u(R c ,z)-u(0,z)\ 2 ) 
(|u(r,L c )-u(r,0)| 2 ) 
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where u(r, z) denotes the deformation field of the vor- 



tex lattice at position (r, z) (with z || B), and 



is the range of the pinning potential.— The transverse 
displacment correlation length 



Rr = 



2$ 



o.A 



1/2 



(3) 



can be obtained, without a priori assumptions, from the 
value of the critical current density. Using the appropri- 
ate parameters (Table [J), one has, for jf y (5K) = 3x 10 9 
Anr 2 , R c = 40 nm in single crystalline PrFeAsOi_j,, and 
R c = 56 nm corresponding to (5K) ~ 1 x 10 9 Am -2 



500 



250 



CO 



-250 - 



-500 - 




-2.0 -1.5 -1.0 -0.5 



1.5 2.0 



FIG. 7: (color online) Hysteresis loops ofthe "self-field", de- 
fined as H a = B_i_//j,o — H a , as measured with a microscopic 
Hall sensor in the center of PrFeAsOi_ H crystal # 7, at the 
indicated temperatures. The values of the onset field H on are 
denoted by the vertical black arrows. 
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FIG. 8: (color online) (a) MOI flux penetration into 
NdFeAsOo.gFo.i crystal # 1, at T = 9.1 K. The three frames 
depict the flux density distribution after the application of an 
applied field of 20 and 50 mT, and the remanent flux after 
removal of the 50 mT applied field, (b) Flux density profiles 
measured across the central part of the crystal. 



in NdFeAsOo.gFo.i. These values are much smaller than 
the intervortex spacing at 300 Oe, at which the data in 
Figs. El El and[9]were obtained. The pinning-induced dis- 
placement of each vortex is thus independent of that of 
neighboring vortices. In this so-called single vortex pin- 
ning limit, one may now estimate the longitudinal dis- 
placement correlation length as 



2 c- 



1/2 



(4) 



one finds L c w 20 nm and 10 nm for PrFeAsOi_j, and 
NdFeAsOo.gFo.i, respectively. This length largely ex- 
ceeds the spacing ofthe FeAs planes, which clearly estab- 
lishes pinning as being in the three-dimensional single- 
vortex (3DSV) limit^ia 

From here on, we show that the critical current den- 
sity in the (1111) iron oxypnictide superconductors can 
be understood as arising from mean-free path varia- 
tions induced by the dopant atoms, oxygen vacancies 
in the case of PrFeAsOi_j,, and F ions in the case of 
NdFeAsOo.gFo.i. The pinning force of a single defect 
is expressed as f p ~ 0.3g(pD)£o (ctr/7r£ 2 ) (£o/£)j where 
a tr = 7rD 2 is the transport scattering cross-section, D v 
is the effective ion radius, and g(po) is the Gor'kov func- 
tion. The disorder parameter po = Tivp /2-kTJ, ~ £o/^ 
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Eq. (10) 
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FIG. 9: (color online) Local values of the critical current den- 
sity j c in the three regions of the NdFeAsOo.gFo.i crystal # 
1, outlined in the central panel of[9ji. Also shown are the val- 
ues of jc V and j'c(O) determined for both investigated crys- 
tals. The first are compared to Eq. ((5} using rid = 1.5 x 10 27 
m~ 3 and D v = 0.9 nm (lower drawn line), while the latter 
are fit to Eq. (|10[1 with respective parameter sets (T c = 37 K, 
m = 6xl0 21 nF 3 , f p , s = 0.1e ) and (T c = 35 K, n, = 2xl0 22 
m~ 3 , /p, s = O.leo). 



with i>f the Fermi velocity, Z the mean free path, and 
£o ~ 1.35£(0) the (temperature-independent) Bardeen- 
Cooper-Schrieffer coherence length i 13 ' The critical cur- 
rent is determined by the fluctuation of the elementary 



pinning force, (/ 2 ), and reads^ 3 . 
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A(0) 
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The numerical factor under the parentheses in Eq. ([5]) 
depends on the precise type of scattering.— Since 
the temperature dependences A(0)/A(T) and £\(T) are 
known from Refs. [2?] and HO (yielding a ~ 2 for 
PrFeAsOi_j, and a ~ 1.5 for NdFeAsOo.gFo.i), one is 
in the position where a full consistency check of both 
the magnitude and the temperature dependence of j c is 
possible^ 

In the case of PrFeAsOi-y, Eq. {5), we start from the 
hypothesis that O vacancies are responsible for the lion's 
share of flux pinning. The ion radius D v = 1.46 x 10 10 
m. Inserting this value into Eq. [5] reproduces the low- 
temperature value jf v = 3 x 10 9 Am~ 2 with the single 
free parameter rid ~ 1.5 x 10 27 m~ 3 . This nicely corre- 
sponds to 0.1 O vacancy per formula unit (half a unit 
cell of volume 65 A 3 ). Eq. {5} reproduces the low-T 
temperature dependence of the critical current density 
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FIG. 10: (color online) PrFeAsOi-y: critical current versus 
doping level, as determined from the phase diagram of Ref.|45l. 
Open circles indicate the critical temperature, filled diamonds 
indicate the critical current density as measured in various 
locations in the four different crystals of Fig. [S] The dashed 
and drawn lines indicate the critical current density expected 
from weak collective pinning by oxygen vacancies, Eq. |(SJ, 
at T = and T — 5 K, respectively (here we suppose that 
A^ 2 oc T c , as in Ref.|47|. 



in the high-j c regions, and the temperature dependence 
over the full range from 5 K to T c in the low-j' c regions. 
The spatially more homogeneous contribution to the crit- 
ical current density of oxygen deficient single crystalline 
PrFeAsO i -j, is therefore well-described by pinning by O- 
vacancies by the 6k mechanism. 

In the case of NdFeAsOo.gFo.i , the analysis is hindered 
by our ignorance of the effective scattering cross-section: 
doping is through chemical substitution, not oxygen de- 
pletion. If one adopts the view that F substitution is 
at the origin of pinning, one has rid ~ 1.5 X 10 27 m~ 3 
for our average doping level. To reproduce the value of 
the measured low-T j c w 1 x 10 9 Am -2 then requires 
<Jtr = 1.5 x 10~ 20 m 2 , corresponding to an effective de- 
fect radius of 0.9 A (this can be compared to the F ion 
radius of 1.3 A). The temperature dependence of jjf v is 
again very well described by Eq. (0. It is not quite as 
strong as in PrFeAsO i_ y , an effect that can be attributed 
to the different T— dependence of the penetration depth 
[A(T) nearly perfectly follows A~ 2 ~ (1 — t 2 )] and of the 
anisotropy ratio (e\ seems to be nearly independent of 
temperature in NdFeAso.gFo.i).— 



B. Spatial variations of j c and link with doping 

Both investigated (1111) compounds show spatial vari- 
ations of both the critical temperature T c and the low- 
temperature critical current density j c . It is tempting to 
correlate the two: knowing the temperature dependence 



8 



of both the superfluid density A~ 2 and the anisotropy 
ratio, as well as the evolution of the respective T c vs. 
doping phase diagrams^^, Eq. (|5|) predicts what the de- 
pendence j c (T c ) should be. In the case of PrFeAsOi_j,, 
our measurements yield sufficient statistics for the ex- 
pected increase of j c with T c to be, indeed, observed. In 
the considered portion of the phase diagram, the more 
vacancies are added, the higher T c , but also the stronger 
the pinning. Fig. [TU] shows a compilation of critical tem- 
peratures and low-temperature critical currents of all in- 
vestigated regions in all our PrFeAsOi_j, crystals. The 
experimental data follow the dependence of the low tem- 
perature j c as this follows from Eq. ((5|) , even though this 
dependence is weak. The contribution to this dependence 
via rid, arising from the addition of oxygen vacancies, is 
actually weaker than the expected contribution from the 
doping dependence of the the superfluid density, which 
we have assumed to follow the relation A~ 2 cx T c %^ Sig- 
nificant scatter due to the strong pinning contribution 
remains in Fig. [TUJ which we shall attribute to the pres- 
ence of doping inhomogeneity on the 10 - 100 nm scale 

In the framework of weak collective pinning, the ob- 
served spatial variation of j c in NdFeAso.gFo.i would, 
if attributed to the macroscopic variation of the dopant 
atom density , correspond to a variation of the doping 
level of x = 0.1 ± 0.03, within a given single crystal. The 
concomittant T c variation would be from 26 K to nearly 
50 K, which is not what is observed by DMO. Moreover, 
and contrary to the observation in PrFeAsOi_j,, the crit- 
ical current density of the investigated crystal is larger in 
areas with low T c , both as far as different regions of crys- 
tal # 1 are concerned, as the observed differences between 
crystals # 1 and 2. In the absence of sufficient statistics, 
we tentatively ascribe this behavior to the presence of 
strong background pinning. 



C. Strong pinning background 
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FIG. 11: (color online) (a) PrFeAsOi- H : Double logarithmic 
plot of the critical current versus magnetic field. The drawn 
lines show the power B~ s ^ 8 expected from the strong pinning 
contribution, (b) ibid , for NdFeAsOo.gFo.i. Drawn lines 
show model fits to Eq. 1101 



As described in section IIII1 the spatial variation of the 
critical current density in single crystalline PrFeAsOi_ y 
is reflected in the temperature dependence of j c , higher 
local j c corresponding to the presence of a break in the 
temperature dependence. Also, the higher local critical 
current densities are responsible for the low-field j c -peak 
observed in Fig. which cannot be explained within the 
single-vortex collective pinning framework. There must 
therefore be supplementary sources of pinning, inhomo- 
geneously distributed throughout the samples, with a 
temperature dependence that is weaker than that of the 
weak collective pinning described above. 

The field dependence of the associated critical cur- 
rent density, a plateau, followed by a power-law decrease 
j c oc _B -5 / 8 , is in very satisfactory agreement with the 
theory of strong pinning developped in Refs. [43II4SL In 
the presence of a density rii of strong pins of size larger 



than the coherence length, one has^ 

1/2 1/2 • / e c \ 3 / 2 



Jc(0) 
fc(B) 



so 



(B < B*) (7) 



1 ■ / f t \ 9 / 4 /rh \ 5 /8 



.5/4.1/2 



Sab 



The 



crossover 



so 



field 



B 



(B > B*) . (8) 



B* 
is 



determined 



0.74^ 2 $ (71.,/uf 5 (f P ,sU/eof 5 
as that above which the so-called vortex trapping area of 
a single pin is limited by intervortex interactions^ The 
identification of the experimental j c (0) with Eq. (JTJ) , and 
of the power-law decrease with Eq. ([5]), allows for the 
determination of the elementary pinning force f p s of a 
strong pin from the ratio [dj s c (B) / 'dB' 5 ^] / [j c s (0)]" 2 . It 
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is found that / P , s (0) = 2 x 10" 13 N for both investigated 
compounds, with a temperature dependence coinciding 
with that of the superfluid density. Hence, we find a 
measured f p>s ~ O.leo- The density of strong pins can 
be straightforwardly estimated from B*: n, « 1 X 10 21 
m~ 3 for PrFeAsOi-y, and rn w 6 x 10 21 m~ 3 and 
« 2 x 10 22 m- 3 for NdFeAsO . 9 F .i crystals # 1 and 2, 
respectively. 

These data can be compared to the results of TEM 
observations. The first candidate strong pins are ex- 
tended (nm-sized) pointlike inclusions or precipitates, 
such as observed in Fig. [T}d. Assuming such defects 
to be non-superconducting, one would have f PtS ~ 
e (A/4£ Q &) In (1 + L> 2 /2£ 2 b ). Typical observed defect 
dimensions are Dj w 2 - 5 nm, yielding / PiS ~ 0.1 - 
l.leo at low temperature. Therefore, the smaller defects 
of radius 2 nm might do the job, were it not that the tem- 
perature dependence expected for such voids is at odds 
with experiment. 

Next, the observed undulations of the FeAs layers im- 
pose an intermittant bending of vortex lines as these 
move through the crystal lattice. The necessary force 
to produce this bending can be estimated as the product 
of the line tension £ 2 £o and the variance (5a) 2 of the tilt 
angle; here, a corresponds to the buckling angle. Such a 
mechanism would yield the experimental temperature de- 
pendence of f Pt g, but, at 10 _4 £o, grossly underestimates 
the measured elementary force. 

Third, the higher strong pinning critical current den- 
sity observed for lower doped NdFeAs(0,F) could be 
linked to the observation of phase coexistence in the un- 
derdoped state of this materials 9 Without going as far 
as invoking the presence of nm-scale magnetically or- 
dered regions in our crystals, the idea of phase coexis- 
tence suggests that there are spatial fluctuations of the 
dopant atom density on the scale of several nm. The en- 
suing dispersion of weakly superconducting regions with 
critical temperature T c — ST C inside a more strongly su- 
perconducting matrix would certainly lead to flux pin- 
ning. Its description would be similar to that of non- 
superconducting precipitates, but with a smaller pinning 
energy, a vortex passing through an area of lower T c gain- 
ing only a fraction ST C /T C of the condensation energy 
£o/4£ 2 . Assuming the condensation energy to be propor- 
tional to the critical temperature, the pinning force can 
be written as 



inserted] and ((SJ), 



Jp,. 



A 

4£at 



In 1 



2& 



(9) 



with t = T/T C) and I = Tj (T c - 5T C ). For small spatial 
variations of the critical temperature, e.g. 8T C /T C ~ 0.05 
or ST C w 1.5 K, and A ~ 5 - 10 nm, Eq. © nicely 
mimics the measured temperature dependence f ps (T) ~ 
£o(T). As shown in Figs. [51 E3 and[TjTb), the total critical 
current density, obtained by summing Eqs. [with © 



3c 



(10) 



is also in good agreement with experimental observations. 
One is thus lead to the conclusion that, in addition to the 
macroscopic inhomogeneity of doping level, there also ex- 
ists an inhomogeneity on the nano-scale, much similar to 
that reported by Yamamoto et al. in Ba(Feo.gCoo.i)2As2. 
However, the doping level modulation, necessarily of the 
order of the T c -variation, ST C /T C ~ 0.05, that explains 
the strong pinning contribution, is far too small to sup- 
port any claims of phase coexistence in the underdoped 
(1111) pnictides investigated here. If similar disorder 
should exist for smaller doping levels, near the super- 
conductivity onset, one would have 5T C ~ T c , and a near 
certain coexistance of magnetic and superconducting re- 
gions. This is a premise that needs further investigation. 

For completeness, one may also contemplate surface 
roughness as a source of flux pinning42r— The critical 
current density is then determined by the force needed 
to push a vortex line out of a surface trough or across a 
ridge, and reads, in the limit of small magnetic field s 43 i 52 



J 



TV 



7T£o 5d 



~ D 2 



(11) 



Here, d would be the crystal thickness, D the spacing be- 
tween surface defects or troughs, and 8d the typical ridge 
height, that is, the variance of the thickness. In Refs.l49i- 
IHlL the ratio 5d/D = sin# c is interpreted as the sine of a 
"contact angle" 8 C . In the Mathieu-Simon model 4 ^— the 
field dependence is expected to correspond to that of the 
vortex chemical potential, i.e. the equilibrium magneti- 
zation. This is not observed. Moreover, if one reinter- 
prets the experimental f p<8 and n, ~ 2/dD 2 in terms of 
surface pinning, one finds a ratio of ridge height to ledge 
width Sd/D ~ 2, for a ledge separation of ~ 20 nm. 
Such a high aspect ratio would mean that the surface de- 
fects are located on the crystal edge, since the alternative, 
cracks on the surface, are not observed. Strong pinning 
by impurities, located in surface regions only, leads to the 
same dependences (|7I8[) . but with 3 x 10 16 defects m -2 . 



D. Fishtail effect and phase diagram 

The knowledge of pinning parameters of the (1111) 
superconductors under study allows one to confront fea- 
tures of the mixed-state (B, T)-phase diagram with theo- 
retical models. In particular, the fishtail effect at B on (T) 
was attributed to a crossover in vortex dynamics as, with 
increasing magnetic field, one leaves the single vortex 



53.54 



or 



pinning regime for the bundle pinning regime 
the occurence of a first order phase transition from an 
ordered, "elastically pinned" low-field vortex phase, the 
so-called Bragg-glass^ to a high field disordered phase 
characterized by the presence of topological defects i 56 i 57 
The latter scenario has been unambiguously verified in 
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FIG. 12: (color online) (B,T) vortex matter phase diagram 
for (1111) iron pnictide superconductors. (Red) Circles in- 
dicate measurements on PrFeAsOi-y, while (blue) squares 
show results on NdFeAsOo.gFo.i- Closed circles show the 
irreversibility field Bi rr (T) measured from the onset of a 
third harmonic reponse from ac Hall probe magnetometry; 
open (blue) squares show the screening onset data of Ref. [29 . 
Peak effect onset fields for both compounds are indicated 
by barred squares (NdFeAsOo.gFo.i) and open (red) circles 
(PrFeAsOi-y). Dotted lines show the single-vortex to bundle 
pinning crossover described by Eq. (|12|l , while dashed-dotted 
lines indicate the order-disorder field described by Eq. (I13p . 



the high temperature superconductors YE^CusC^-a 
and I^Si^CaC^Os+^r^— > m the cubic superconduc- 
tor (Ba,K)Bi0 3 ^ in NbSe^ as well as in MgB 2 ^i 

In the first case, the onset field B on should coincide 
with the single-vortex- to bundle pinning crossover field 
Bsv, determined by the equality of R c [see Eq. @] and 
the vortex spacing ao: 



B S \ 



405, 
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(12) 



Inserting the experimentally obtained jf v into Eq. (JTJ 
yields the dotted lines in Fig. [T2J Clearly, while 
the experimental B on data for more strongly pinning 
PrFeAsOi_j, lie below those for more weakly pinning 
NdFeAsOo.gFo.i, Eq. [T2l predicts otherwise. Therefore, 
even if the peak effect onset lies in the vicinity of the 
single- vortex to bundle pinning crossover, it cannot be 
directly associated with it. 

On the other hand, the vortex ensemble can undergo 
a structural transition whereby it lowers its energy by 
adapting itself more efficiently to the underlying pinning 
potentiel, at the expense of the generation of topological 
defects^ - — In the absence of a theory for this order- 
disorder transition of the vortex lattice, a Lindemann- 



like criterion was developped in Refs. I55ll59l and [6^ in 
order to, at least, estimate its position in the (B,T)- 
plane. The Lindemann approach considers that topolog- 
ical defects can be generated when pinning is sufficiently 
strong to provoke the wandering of vortex lines outside 
their confining cage formed by the nearest neighbor flux 
lines. The different result o 59 i 60 have been summarized 
in Ref. l57l In the regime of single vortex pinning, rele- 
vant for collective pinning in the (1111) compounds, the 
position of the order-disorder transition is given by 
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AU SV°OD 
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where b on = B d/B c2 , b S v = B S v/B c2 , c L ~ 0.1 is the 
Lindeman number, A is a numerical constant, t = t/T c , 

and F T {t) = 2t (Gi/1 - t 2 ) 1/2 . The use of the param- 
eters of Table HI the experimentally measured jf v , and 
A = 4 yields the dashed lines in Fig. [T2] These show 
more than satisfactory agreement with the experimen- 
tally measured positions of B on . We therefore conclude 
that, most likely, a bulk order-disorder transition of the 
vortex lattice lies at the origin of the peak effect in (1111) 
pnictide superconductors. However, more work, espe- 
cially on vortex dynamics and possible hysteresis associ- 
ated with the transition, should be performed to ascer- 
tain this. 



V. CONCLUSION 

It is found that superconducting iron pnictide single 
crystals show significant spatial variations of both the 
critical temperature T c and the critical current density 
j c . Variations of these quantities on the macroscopic 
scale, from several to several hundred /im, are at the 
origin of a smearing of globally measured properties, and 
notably of the width of the superconducting transition. 
This implies the necessity of local measurements, such 
as magneto-optical imaging or Hall-probe magnetometry, 
to extract superconducting parameters. From such local 
measurements, it is found that the critical current in iron 
oxypnictide superconductors of the (1111) family of com- 
pounds arises from two distinct contributions. The first 
is weak collective pinning by dopant atoms or vacancies, 
vortex lines being pinned by the small scale fluctuations 
of the local dopant atom density. The pinning mecha- 
nisme is identified as being due to mean-free path vari- 
ations in the vortex core (5k mechanism). This means 
that dopant atoms should also be effective quasi-particle 
scatterers. The second pinning contribution manifests it- 
self at low fields. The corresponding critical current con- 
tribution can be completely parametrized by the strong 
pinning theory developped in Refs. [43ll4sL which means 
that extended defects are at its origin. An analysis of the 
magnitude and field-dependence of this strong pinning 
contribution shows that spatial variations of the doping 
level on the scale of several dozen to one hundred nm may 
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be at stake. These variations do not support the possi- 
ble coexistence of the anti- ferromagnetic metallic and the 
superconducitng phases. Finally, we contend that a bulk 
order-disorder transition of the vortex ensemble is at the 
origin of the "fishtail" or peak effect observed in the crit- 
ical current in sub-T fields. 
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